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Abstract 

A theory for the prediction of the size dependence of torsional rigidities of nanosized 
structural elements is developed. It is shown that, to a very good approximation, the 
torsional rigidity (D) of a nanosized bar differs from the prediction of standard continuum 
mechanics (Dc) as {D — Dc) / Dc = Aho/a where ^ is a non-dimensional constant, a is the 
size scale of the cross-section of the bar and ho is a material length equal to the ratio of the 
surface elastic constant to the bulk elastic constant. The theory developed is compared 
with direct atomistic calculations ( "numerical experiment" ) of the torsional rigidity bars 
made of several FCC metals modeled using the embedded atom method. Very good 
agreement is obtained between theory and simulation. The framework presented here can 
aid the development of design methodologies for nanoscale structural elements without 
the need for full scale atomistic simulations. 



1 Introduction 

The demand for smaller and faster devices have encouraged technological advances resulting 
in the ability to manipulate matter at micro- and nanoscales that have enabled the fabrication 
of micro/nanoscale electromechanical systems. While MEMS technology is now a well estab- 
lished area, nanoelectromechanical systems (NEMS) have recently made their appearance in 
literature (see, for example, Roukes (2000)). A key feature of NEMS that make them attrac- 
tive are their high fundamental frequencies while affording small force constants. Nanosized 
bars and tubes are produced in a variety of ways with materials such as SiC, M0O3 and C 
(carbon nanotubes) (see Yakobson and Smalley (1997), Terrones et al. (1999), Sheehan and 
Lieber (1996)). These nanosized elements have found many technological uses, for example, as 
probes of scanning probe microscopes (Dai et al. 1996), in altering properties of bulk materials 
in the form of whisker additions (Kuzumake et al. 1998) etc. 

In the recent past, several groups have reported studies on the mechanical behaviour 
nanosized bars and nanotubes. Wong et al. (1997) performed experiments on SiC beams while 
Poncharal et al. (1999), and more recently Gao et al. (2000) have reported experiments where 
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the elastic modulus of carbon nanotubes is measured using dynamic techniques. A careful 
study of these reports show that the elastic moduli of such nanosized structural elements 
depend on their size. Attempts to explain the size dependent behaviour has been through 
direct atomistic computer simulation of these structures (Robertson et al. 1992, Garg and 
Sinnott 1998, Garg et al. 1998). These studies have reported size-dependent elastic moduli 
computed from atomistic simulations. 

Miller and Shenoy (2000) developed a simple model to explain the size dependence of the 
elastic rigidities of nanosized structural elements. These size dependences were attributed to 
the heterogeneities in the atomic environments introduced by the bounding free surfaces of 
the structural elements. Thus, as the size of the structure becomes smaller the presence of 
surfaces have to be accounted for in modeling strategies. Developing this premise. Miller and 
Shenoy (2000) showed that the differences between the rigidities (D) of these small elements 
and those predicted by continuum mechanics (Dc) can be expressed as 



where A is a nondimensional constant that depends on the geometry of the structure, a is the 
size scale of the structural element (for example, the cross-sectional width of a bar) , and ho is 
a material length that is the ratio of the surface elastic constant of the bounding surfaces of 
the structure and the bulk elastic constant of the material. Thus, the size dependence of the 
rigidities can be predicted by obtaining the material parameter ho and the nondimensional 
constant A. Typically, ho can be obtained from a small atomistic simulation and A can be 
calculated analytically. Thus the need for full scale atomistic simulations of structures (which 
is an expensive proposition) is obviated. Miller and Shenoy (2000) apphed this model to study 
the elastic properties of nanosized bars, plates and beams and demonstrated the strength of 
the model by comparison with direct atomistic calculations. 

An important mode of deformation of bar-like structures is torsion; for example, the probe 
of a scanning probe microscope is subjected to bending and torsion. It is therefore important 
to develop a model for the torsional rigidities of nanosized elements; this is the aim of the 
present paper. An augmented continuum theory of torsion accounting for the presence of 
free surfaces is developed and the size dependence of the rigidity is derived analytically. A 
perturbative scheme is developed for solving the resulting boundary value problem which 
provides a simple method to evaluate the size dependence of the torsional rigidity. It is shown 
that, in general, the size dependence of the rigidity D is of the form 



where A, B, etc. are non-dimensional constants and = S/G where S is the surface shear 
modulus and G is the bulk shear modulus. The perturbative method also provides a general 




(2) 



framework for the calculation of the nondimensional constants A, B, . . . that depend only on 
the geometry of the cross section of the bar. For the case of square bars of side 2a, the theory 
provides that size dependence of torsional rigidity to be very well approximated by 

^ 4^. (3) 
Dc a 

These theoretical results are then compared with direct atomistic simulations (which serve as 
numerical experiments) of torsion of square bars of FCC metals. The methodology required 
for the atomistic simulation of torsion is also developed here. The agreement theory and 
numerical experiment (atomistic simulations) is excellent. 

The paper is organised as follows. The next section contains the augmented theory of 
torsion. Section 3 contains the details of simulation methodologies, the results of which are 
reported and discussed in section 4. The paper in concluded in section 5 in which several 
important directions of future work are identified. 



2 Theory 

2.1 Augmented Continuum Theory 

It was shown by Miller and Shenoy (2000) that the elastic properties of nanoscale structural 
elements such as plates, bars and beams can be explained using a augmented continuum theory 
that accounts for the energetics of deforming inhomogeneities such as surfaces and corners, the 
effects of which are significant when the length scale of the structure approaches the atomic 
scale. Several authors had previously utilized continuum theories of solids with surface effects 
(Gurtin and Murdoch 1975, Rice and Chuang 1981, Cammarata 1994) to study a variety 
of problems ranging from diffusive cavity growth in stressed solids to stability of stressed 
epitaxial films. The formulation outlined in Miller and Shenoy (2000) is briefiy recapitulated 
here for the sake of completeness and to set the notation. 

The body B, described by coordinates Xi, considered in the augmented continuum theory 
is bounded by a surface §. It is assumed that the surface S is piecewise fiat (this assumption 
eliminates the need to consider contravariant and covariant components of surface tensors) 
and is described by coordinates for each fiat face. The bulk stress tensor in the body S is 
denoted by aij and the surface stress tensor by r^fl. Mechanical equilibrium of a bulk material 
element implies that the bulk stress tensor satisfies (with no body forces) 

aijj = 0. (4) 

Equilibrium of a surface element necessitates that 

ra/3,l3 + fa = (5) 



(6) 



where rii is the outward normal to the surface, fa is the negative of the tangential component 
of the traction ti = (JijUj along the a direction of surface S, and is the surface curvature 
tensor. The assumption of the piecewise flat surfaces implies that the surface curvature 
vanishes everywhere along the surface except at corners and edges which have to be treated 
separately. It must be noted that the assumption of piecewise flat surface is merely for the 
sake of mathematical simplicity; the present theoretical framework is valid for curved surfaces 
as well. 

The kinematics of the body is described by the displacement fleld Ui deflned at every point 
in the body. The strain tensor eij in the body is obtained using a small strain formulation as 

^ij = \ {uid + Uj^i) ■ (7) 

The surface strain tensor is derived from the bulk strain tensor e^j such that every material 
fibre on the surface has the same deformation whether it is treated as a part of the surface or 
as a part of the bulk, i. e., the surface strain tensor is compatible with the bulk strain tensor. 

The final ingredient of the augmented continuum theory is the constitutive relations that 
relate the stresses to strains. The bulk is considered to be an anisotropic linear hyperelastic 
solid with a free energy density W defined as 

^{^ij) = -^CijkieijSki (8) 

and the stresses are derived as 

dW 

o'ij — = Cijki€ki (9) 

where Cijki is the bulk elastic modulus tensor. In this framework the bulk free energy vanishes 
for the unstrained solid. Constitutive relations for the surface stress tensor are more involved. 
The surface stress tensor is related to the surface energy 7 as 

Ta/3 = 7<^a/3 + ^ , (10) 

a relation which is generally attributed to Gibbs (Cammarata 1994). The surface stress tensor 
can be expressed as a linear function of the strain tensor as 

Tal3 — Tal3 + Sal3^S<^"/5 (H) 

where is the surface stress tensor when the bulk is unstrained (obtained from (10) with 
^a(3 = 0) and Sap-ys is the surface elastic modulus tensor. This is an important quantity in 
that the size dependence of elastic properties will be shown to be determined by the ratio of 
a surface elastic constant and the bulk elastic constant. The constitutive constants Cijki and 
Sai3js s-re external to the augmented continuum theory; in this paper these are determined 
from atomistic models of the materials considered. 




Figure 1: Cross section of bar considered in the augmented continuum theory of torsion. 



2.2 Augmented Theory of Torsion of Bars 

The continuum theory of torsion of bars is a part of classical theory of elasticity attributed 
to St. Venant and Prandtl and is treated in much detail by Sokolnikoff (1956). In this section 
we develop a theory of torsion based on the augmented theory of last section which includes 
surface effects. Corner effects are neglected in this treatment. 

Fig. 1 shows the cross-section of the bar S bounded by a surface S. Attention is restricted 
to simply connected cross- sect ions for the sake of simplicity. The outward normal to § is 
denoted by n and s is the tangent vector to § with n ■ s = 0. Application of a torque T to 
the bar produces a twist per unit length a; the aim of the analysis is to obtain a relationship 
between T and a. The kinematics is described by the displacement field as 



Ml(xi,X2,X3) = a(f){x2,X3) 
U2{Xi,X2,X3) = -aXiXs 
U3{Xl,X2,X3) = aXiX2 



(12) 



where is the warping function. The only nonvanishing strain components derived from these 
displacements are 



a f d(f) 

ei2(X2,X3) = - ( — Xs 



The bulk material is assumed to be a linear elastic solid and thus the bulk stresses are related 
to the strain via 



(Ti2 = 2Gei2, (Ti3 = 2Ge 
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(14) 



where G is an appropriate shear modulus (for example, if the '1' direction corresponds to the 
(100) direction in an FCC crystal, G = C44). The only nontrivial equilibrium equation in the 
bulk is 



dai2 d(Ti-i 

— — H — = 0. 

dx2 dxs 



(15) 



Eqn. (15) is identically satisfied on the introduction of a stress function ip (after Prandtl) such 
that 

(7i2 = Ga——, (713 = -Ga——. (16) 

On substituting the expressions for the stresses (16) in (14) and using (13) it is found that 
the stress function satisfies 

dxl dxl 

in the bulk. 

The surface is parametrised by the coordinate Xi and distance s measured along the tangent 
vector s. The pertinent component of surface stress is where s stands for the s direction 
in fig. 1 and satisfies the surface equilibrium equation (5) 

^ + A^0. (18) 

The second surface equilibrium condition (6) is identically satisfied since Kq/? — 0. The "surface 
body force" /i is given as 

fi = -ti = -(cTiana + ai^ris) 

^ f dil> dip 
= -Ga ^— na - ^— ^3 
\ 0x3 0x2 

f dip dip 
= -Ga —S3 + — S2 

\OX3 0x2 

The surface stress is related to the surface strain as 

Tls = 2^61, (20) 

where 5" is an appropriate surface shear modulus (r^^ can be taken to be zero with out loss of 
generality), i. e., 

Tls = 25615 = 25'(-ei2n3 + ei3n2) 
= 25(--n3 + -n2) 
[dip 

= -Sa ^—n2 + -K—n^i 

\OX2 OXs 

Substituting (19) and (21) in (18) the boundary condition on S for -0 is obtained as 



where C is a constant along the boundary. It is assumed that S does not depend on s, i.e., 
all the bounding planes are assumed to be crystallographically equivalent. 

The torsional rigidity of the bar can be computed as follows. The torque T has contribu- 
tions from the bulk "B and the surface S and is given as 

T = (X2cri3 - X3(7i2)dX2dX3 + / Tis{x2n2 + X3n3)ds 
J2 J$ 

= 2Ga J ^dx2dx3 -Ga + ^1^) ^^^^ 

Evidently, the constant C can be chosen to be zero. The torsional rigidity D can be obtained 
by solving the mixed boundary value problem 



Gdn 

and substituting for %Ij in the expression 

D^2G 



dxl dxl 



1/; + ——- = on 



(24) 



/ Vdxadxa. (25) 



Clearly, the theory reduces to the standard theory of torsion when S is set to zero. It is 
also clear that the key parameter that determines the atomistic surface effects that affect 
the torsional rigidity is the material length- scale ho defined by the ratio of the surface shear 
modulus S and the bulk shear modulus G, i. e., by 

^0 = |. (26) 

In the next section, a perturbative solution to the boundary value problem (24) will be devel- 
oped and a general formula for the size-dependent torsional rigidity will be obtained. 

2.3 Perturbative Solution 

The general perturbative solution will be developed in a non-dimensional form. To this end, 
the geometry of the cross section of the bar is assumed to be characterised by a length scale 
a (for example, if the cross section is a square, a can be chosen as one half of the side of the 
square). The following non-dimensional quantities are introduced 

* = 4 e=£?, , = £2. (27) 

In terms of these nondimensional quantities, the boundary value problem (24) can be recast 
as 

V^* = -2 in S 

(28) 

a* 

* + = on § 

on 



where = ^ + ^ and the rigidity is 

D 

with 



2 / ^dCdr/, (29) 



= 2 / ^odCdr? (31) 



/3 = ^, (30) 
a 

the nondimensional parameter that governs the extent of surface effects. When (3 — the 
boundary value problem is solved by the nondimensional stress function ^-nd the torsional 
rigidity is 

Dc 

where the subscript c is used to denote that standard continuum value of the torsional rigidity. 

When /5 7^ 0, the nondimensional stress function can be expanded in a perturbative series 
in /3 as 

* = ^'o + /3*i + /9'*2 + --- (32) 
Since ^ solves the boundary value problem (28), it follows that 

/3V^*i + /3^V^*2 + ... = (33) 

in S and 



on §. Since the perturbative expansion (32) is valid for all values of (33) and (34) imply 
that 

V^^i = on S with *i = — — ^ on § (35) 

on 

V^*2 = on S with *2 = — ^ on § (36) 

on 



V^^fc = on S with = ^ on § (37) 

On 

which is a sequence of Dirichlet boundary value problems. Thus, the solution of the boundary 
value problem (28) can be obtained to any desired accuracy in the parameter /3, and the 
torsional rigidity obtained as 

-^^Ao + A,/3 + A2P^ + ... (38) 



where 

Ak^2 f *fcded77 (39) 
is 

are constants that depend only on the shape of the cross-section. Since Dc/Ga^ — Aq, a, 
general expression for the size dependence of the torsional rigidity can be derived as 

For cross-sections of a given shape, the constants Ak can be calculated once and for all and 
the formula (40) along with (26) and (30) can be used to predict the size dependence of the 
torsional rigidity. Atomistic inputs are required only in providing values for S and G. It will 
be seen in the following sections that the value of (3 in real systems is less that 0.5 and the 
most important contribution is from the first term on the right hand side of (40), i. e., 

£^^^p = ^h (41) 

3 Atomistic Simulations 

Atomistic simulations are carried out to validate the theory developed in the previous sections. 
The bars selected for study have a square cross section (fig. 2) made of selected FCC metals 
(Al, Ag, Cu, Ni) such that the 'I'-direction corresponds to the [100] crystallographic direction 
and '2' and '3' directions correspond to [010] and [001] crystallographic directions respectively, 
i. e., the bounding free surfaces are planes of the {100} family. 

Two sets of simulations are carried out. In the first set, the constitutive constants G — C44 
(bulk shear modulus) and S the surface elastic constant for surface shear of the {100} surface 
are evaluated atomistically; these are the parameters required as input to the theory developed. 
In the second set of simulations the torsional rigidity is directly calculated using atomistic 
simulation of nanoscale torsion. These results are then compared with the theoretical results 
of torsional rigidity in the next section. 

The atomistic model used in the present study is the embedded atom method (EAM) 
developed by Daw and Baskes (1984). The elements Ag, Cu, and Ni are modeled using the 
EAM potentials of Oh and Johnson (1988) and Al is modeled with potentials developed by 
Ercolessi and Adams (1994). 

3.1 Determination of Surface Elastic Constants 



The surface elastic constants arc determined as follows. A block of atoms are stacked in an 
FCC crystal lattice such that the coordinate planes are of the {100} family. Periodic boundary 
conditions are imposed in the '1' and '2' directions alone simulating bounding free surfaces of 



Material 


ao (A) 




S (eV/A2) 


ho (A) 


Al 


4.032 


0.229 


0.481 


2.099 


Ag 


4.090 


0.292 


0.248 


0.849 


Cu 


3.615 


0.474 


0.429 


0.906 


Ni 


3.520 


0.808 


0.763 


0.945 



Table 1: Properties of materials calculated using EAM potentials. Oo-lattice parameter, C44- 
bulk shear modulus, ^-surface shear modulus and Hq = S/C44. 

{100} family in the '3' direction. The positions of the atoms are changed to correspond to a 
deformation gradient tensor F = I + eci (E) 62 that produces a simple shear by amount e (/ 
is the identity tensor, and Cj are the orthonormal basis vectors). The potential energy of this 
atomistic system is minimised and the minimised total energy is calculated. To compute the 
surface energy, the elastic energy stored in the bulk is subtracted from the calculated total 
energy. On performing this simulation for various values of e, the surface energy 7 is obtained 
as a function of e, and the surface elastic constant is calculated by numerical differentiation 
of this function. It is noted that to compute the bulk elastic constant G{= C44) a similar 
procedure can be adopted with periodicity imposed in the '3' direction as well. The results of 
these simulations are shown in table 1. 

3.2 Atomistic Simulations of Torsional Response 

Atomistic simulations used to calculate the torsional rigidities are performed as follows. The 
simulation box {0 < Xi < i,\x2\ < a,\x3\ < a) consists of a collection of atoms as shown in 
fig. 2 with periodic boundary conditions imposed in the '1' direction with a periodic length of 
i. The atoms are then displaced according to the rotation tensor 

R{xi) = / + a (63 (g) 62 - €2 ® 63) (42) 

where a is the twist per unit length {a is chosen such that ai <^ 1, so that the infinitesimal 
rotation tensor is sufficiently accurate). With this initial condition, the atomistic potential 
energy of the system is minimised. The energy and force computation during the minimisation 
process is carried forth as follows. To compute the energy of an atom (say i) near xi — i the 
position vector of its periodic neighbour (say j), found near xi = in the box, is required. 
The position vector of the periodic image of atom j is obtained by applying the rotation tensor 
R{t) to the vector x{j) + iei where x{j) is the position vector of atom j in the simulation 
box. Similarly the periodic neighbour j of an atom i near xi — can be obtained by applying 
the inverse of the rotation tensor R{£). This procedure ensures that the periodicity in the '1' 
direction is maintained, while keeping the bar in a twisted position with twist per unit length 
a. 
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Figure 2: Atomistic system used for the simulation of torsional response. 

Using the procedure above, simulations are carried out for various values of a and the 
minimised atomistic total energy E in the simulation box is obtained as a function of a. The 
torsional rigidity is calculated using the formula 



(43) 



The correctness of this procedure is ascertained by choosing various values of i for the periodic 
distance and computing D using the above procedure. It is found that D is insensitive to the 
choice of i. 



4 Results and Discussion 

4.1 Theoretical Results for Square Bars 
4.1.1 Exact Solution 



The exact solution of the boundary value problem (24) for the case of a square bar is obtained 
in this section. The square section is assumed to be of side 2a and the solution for the boundary 
value problem is obtained in the nondimensional form (28) where the square occupies the 



region \^\ < 1, |?7| < 1. A straightforward analysis gives that 

^) = 8 E kl{2kn + 2n2kn) ~ coshA;„ + /3£linhA;J ^^^^ 

where A;„ is the n-th root of the equation 

cos k — Pksink = 0. (45) 
The nondimensional theoretical torsional rigidity (29) is obtained as 

J^ = 64y (l V (46) 

Ga^ ^ kl(2kn + sin 2kn) \ A;„(cosh kn + Pkn sinh kn) J 

For the purpose of comparison with atomistic simulations, the nondimensional warping 
function $ = is also derived: 

oo . ^ 

m ^) = ^ ^ - 8 E fc^(2fc. + sin 2fc„) (cosh + sinh M ^^"^ ^''^ 

4.1.2 Perturbative Solution 

The function is obtained by setting = in the exact solution (44). The function is 
obtained by solving the Dirichlet problem (35) as 

n=oo . , , 

V ^ smg„tanhg„ 

Wi(4, T]) = A (cos Qnt. cosh g„r] + cosh Qnt. COS qnT]) (48) 

„=0 cost! Qn 

(2 7/ — l)7r 

where Qn = . The constant Ai can be calculated as 

\ ^ sill (I., taiiii q„ 
= 64 5^ (49) 

n=l ^" 

The numerical values of Aq (which is calculated from (46) with /3 = 0) and Ai are 

Ao = 2.2492 Ai = 8.9969 (50) 

Thus, 

DKiAo + Ai(5 (51) 

and 

(62) 

Dc Aq a 
A comparison between the exact result (46) and the perturbative result (51) is shown in 
fig. 3. It is evident that the perturbative result is a very good approximation of the exact 
result; even when (3 — 1 (which is much larger than that which would appear in systems 
considered here), the error is only 6.5%. Thus the perturbative result in the form of (52) is 
used in comparisons with atomistic results. 




Figure 3: Comparison of the exact (46) and perturbative (51) results for the torsional rigidity. 

4.2 Comparison of Atomistic and Theoretical Results 

Fig. 4 shows a comparison between the nondimensional torsional rigidity computed atomisti- 
cally and that predicted by the augmented continuum theory (52). Several points may be 
noted: 

1. It is clear that the atomistically calculated torsional rigidity differs significantly from 
the predictions of standard continuum theory. In fact, for an Al bar of width 5nm the 
torsional rigidity is about 50% larger. 

2. The nondimensional difference in torsional rigidity computed atomistically scales very 
closely as (1/a) in all the metals. 

3. The atomistically computed values arc accurately predicted by the theoretical values for 
Al and Ni (within 10%), while the agreement is fair for Ag and Cu (within 30%). 

4. In all cases, the atomistically computed values are larger than the theoretical values, 
although to varying degrees in different metals. 

It may be argued that the reason for the atomistic values being greater than the theoretical 
values is the neglect of corner effects in the theoretical analysis. A simple dimensional analysis 
indicates that corner effects must scale as 1/a^; but the atomistic results scale very closely 
as 1/a. Thus it is clear that corners play a secondary role in the systems considered here. A 
more plausible reason for the difference in the theoretical values and simulation results is the 
assumption that the surface energy 7 depends only on surface strain. In reality the surface 




Figure 4: Comparison of the results of atomistic simulations with theoretical results 
lattice parameter. 



Oo is the 



energy can also depend on the surface curvature strain bafs, i. e., the difference between 
the deformed curvature and the original curvature of the surface. Thus the surface energy 
must be a function of both the surface strain tensor and the surface curvature strain tensor; 
mathematically, 7 = 7(60,^,60^). This idea can be expressed more physically; in the present 
model, the surface is treated as a membrane, while a more physically realistic model will be 
that of a shell with bending stiffness. Additional evidence in support of this argument is that 
in the case of bending of plates treated by Miller and Shenoy (2000) the augmented continuum 
theory differs by the order of 30% with the atomistic simulations. The plates treated did not 
have corners and this effect can only be attributed to the neglect of the bending energy of the 
surfaces. Clearly this problem requires a more elaborate theoretical framework and will be 
taken up for study in future. 

The atomistic simulations of torsion not only allow for the computation of the torsional 
rigidity but also provide data for the warping of atomic planes. From the positions of atoms at 
the configuration of minimum energy, the values of nondimensional warping $ are calculated. 
The atomistically computed values of the nondimensional warping must be independent of 
the value of a (twist per unit length) according to the theory. This is indeed found in the 
simulations, and provides a further check for the theory. The atomistically simulated warping 
displacements are then compared with the theoretically predicted values (47). The result of 
such an exercise is plotted in fig. 5. It is evident that the predicted warping is in excellent 
agreement with the atomistic result. 

5 Conclusions 

A general framework for the prediction of rigidities of nanoscale structural elements has been 
developed and applied to the case of the nanoscale bars in torsion. The key premise of this 
theory is that the heterogeneities present in such systems can be modeled as surface effects in 
an augmented continuum theory. A development of this idea reveals that the material length 
scale that governs the size dependence of the rigidity is the ratio of the surface elastic constant 
to the bulk elastic constant. The augmented theory of torsion developed here is compared 
with direct atomistic calculations of bars of various metals and is found to be satisfactory. 
The author is not aware of any work that reports experimentally measured torsional rigidities 
of nanosized bars - indeed, atomistic simulations are used as numerical experiments. Given 
the advances in nanotechnology, such experiments are expected to be performed in the near 
future. 

The use of this theory is envisaged as follows. The bulk elastic constants and the surface 
elastic constants (for various surfaces) of materials of interest can be calculated and tabulated. 
The expressions for the constants that appear in the perturbative expansion in size dependence 
of the torsional rigidity can be worked out for a host of cross sectional shapes once and for all. 
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Figure 5: Comparison of atomistically simulated and theoretical warping function $. The 
solid lines are the contours of the atomistic result while the dashed lines correspond to the 
theoretical calculation (47). The atomistic result is for Al with a/ao = 5 and the theoretical 
result corresponds to /3 = 0.1. 



A collection of such information will be useful for the designers of nanomechanical systems in 
that the need for direct atomistic simulations of nanoscale structures is obviated. The work 
of Miller and Shenoy (2000) along with the present work provide a complete framework for 
the prediction of rigidities of nanoscale structural elements in extension, flexure and torsion. 

Several points for future work are noted. The atomistic model used in this study, EAM, 
is known to be inaccurate with applied to interfacial properties. This, of course, does not 
invalidate the present work since the parameters used in the theory and the simulation results 
are obtained from the same EAM model. To obtain accurate values of surface elastic con- 
stants more sophisticated atomistic models such as density functional theory (DFT) may be 
applied. Another important point to be noted is that thermal effects are not accounted in the 
calculation of the surface elastic constants. Also, the model developed here treats the surface 
as a membrane while a more realistic model will have to additionally account for the bending 
stiffness of the surface. Results of investigations along the afore mentioned directions will be 
reported in future publications. 
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